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Abstract 

A canonical particle definition via the diagonalisation of the Hamiltonian for a quan- 
tum field theory in specific curved space-times is presented. Within the provided approach 
radial ingoing or outgoing Minkowski particles do not exist. An application of this for- 
malism to the Rindler metric recovers the well-known Unruh effect. For the situation 
of a black hole the Hamiltonian splits up into two independent parts accounting for the 
interior and the exterior domain, respectively. It turns out that a reasonable particle 
definition may be accomplished for the outside region only. The Hamiltonian of the field 
inside the black hole is unbounded from above and below and hence possesses no ground 
state. The corresponding equation of motion displays a linear global instability. Possible 
consequences of this instability are discussed and its relations to the sonic analogues of 
black holes are addressed. 
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1 Introduction 



At present we know two very fundamental and effectual theories in order to describe nature, 
quantum theory and general relativity. But a satisfactory unification of these distinct theories 
is still missing. One possibility to achieve some progress towards this aim is expected to be 
provided by the consideration of quantised fields in given classical space-times. Within this 
treatment the metric plays the role of an external background field. Various investigations 
have been devoted to this topic during the last decades, to mention only some of the most 
important initial papers in chronological order: In 1972 Fulling noticed the non-uniqueness 
of the particle interpretation in curved space-times (which may be regarded as the basis for 
several effects) . Two years later Hawking || found out that black holes are not completely black 
but possess a thermal behaviour caused by the conversion of the initial vacuum fluctuations. A 
short time after this striking discovery Davies H] showed that also a uniformly accelerated mirror 
- treated as a mirror at rest in the Rindler metric - creates a thermal spectrum. In 1976 Unruh 
0| recognised the fact that even a uniformly accelerated observer in the Minkowski vacuum feels 
environed by a thermal bath. Many examinations have been accomplished since these basic 
papers, see e.g. [Q— [3T| and references therein. Now black holes are very interesting touchstones 
in order to test candidates for the theory that unifies general relativistic and quantum aspects. 
The representations of black holes within the underlying theory are expected to reproduce their 
main properties. 

The main intention of this article is to provide a canonical approach to quantum field theory 
in specific curved space-times and a related particle definition, together with an investigation 
of the consequences of this formalism. 

This paper is organised as follows: In Section |] we set up the equations which describe the 
quantum field under consideration together with the assumptions necessary for the particle 
definition. Based on methods of functional analysis we propose in Sec. ^| a canonical approach 
to the particle definition via the diagonalisation of the Hamiltonian. This procedure is applied 
to some flat space-times in Section ^ in order to elucidate the underlying mechanism. Sec. |5] 
considers the space-time of a black hole and Section |5.1| presents the particle definition for the 



exterior domain. In Sec. |5.2| we examine the quantised field inside the black hole and deduce the 
unstable behaviour of its time-evolution. Possible consequences of this instability are discussed 
in Section |5.3| and its implications for the sonic analogues of black holes are pointed out in 



Section |5.4| . We shall close with a summary, some discussions and an outline. 
Throughout this article natural units with G = h = c = k^ = l will be used. Lowercase 
Greek indices such as /i, v vary from (time) to 3 (space) and describe space-time components. 
Lowercase Roman indices i,j range from 1 to 3 and label only the space. (For both we employ 
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the Einstein sum convention.) Uppercase Roman indices /,J 6 N may assume all natural 
numbers while uppercase Greek indices like T, A may be more general, e.g. continuous. 



2 Equations of motion 

We consider a minimally coupled, massless and neutral (i.e. real) scalar field $ whose propa- 
gation in the space-time (M, g^ v ) is described by the action 



a= l d L z*-^-<r(d^)(a v *), (i) 



with g = det(g fJiU ). Possible potential terms like a mass term m <& or a conformal coupling 
term TZQ 2 /6 (where TZ = TVt denotes the Ricci scalar) do not alter the main conclusions, see 



Sees. |3.2| and |5.2j . The same holds true for a charged and thus complex field $ and $*. For 
reasons of simplicity and considering the scalar field $ as a model for the photon field we restrict 
ourselves to the most simple action in Eq. (0). 

Provided that the spatial surface terms arising from the integration by parts vanish the variation 
of the action 5 A = leads to the Klein-Fock-Gordon equation 

□$ = d, grd v Q) = . (2) 

The corresponding inner product is defined by 

=i J dXr^dpt (3) 
s 

with ip d[i 4> = d^cp — (pd^ip. In this definition the surface element effi^ already contains 
factors like y/—g-£ with the result dH^dx^ = d 4 x yj—g. For functions which fulfil the Klein- 
Fock-Gordon equation Dip = Dtp = the inner product is independent of the special surface S, 



cf. for instance [29|. To see that, one has to use Gauss' law 



dS M A 1 = J <Tx V M A" (4) 
and again to require vanishing spatial surface terms. 
2.1 Preconditions 

Now we have to specify the assumptions which are necessary for an appropriate particle defini- 
tion. At first we demand a strongly causal space-time M. This condition ensures the essential 
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physical principle of distinguishing cause and effect and forbids (for instance) the occurrence 



of closed time-like curves, cf. [33 



As another requirement we impose a static metric of the space-time M 

ds 2 = g^dx^dx" = g 00 (r)dt 2 + g ij (r)dx t dx j . (5) 

These two assumptions allow to factorise the space-time M = M Cg> G into time t EM. and space 
r E G with an open domain GCK 3 . The Killing vector corresponding to the time translation 
symmetry permits the definition of a conserved energy. This fact is substantial for a physical 



reasonable particle definition, see also Sec. E72] below. On the other hand, the selection of 



a particular Killing vector refers to a class of associated observers whose time evolution is 
generated by this vector field. In general different Killing vectors generate distinct particle 
definitions applying for the different observers, see Sec. ^]3| below. 

The third precondition we need is called non- degenerated signature. This fixes the signature of 
the metric inside the domain G 

VreG : g 00 (r)>0 ; (^(r)) y . < . (6) 

In the latter inequality (gij(v))- has to be understood as a matrix (and not as the single 
components), i.e. Vp E M 3 : p 2 > — > p%g^Pj < 0. Both quantities (g o an d gij) are continuous 
and regular inside G but may diverge or vanish by approaching the boundary dG. This might 
be the case for a horizon situated at dG. (A g 00 -component of the metric that vanishes over a 
finite volume would create primary constraints, see e.g. f32fl .) 



As demonstrated above the possibility of performing the integration by parts is a really im- 
portant issue. Accordingly, our last assumption is a physical complete region G. This simply 
enforces the vanishing of the surface terms. The occurrence of such boundary contributions 
always indicates the interaction with a system behind the surface. Such a region would not be 
physical complete. It should be noted here that the validity of the integration by parts also 
includes periodicity in angular coordinates, such as f(ip) = /(<£> + 2tt). For the above specified 
space-time M = M <S> G the spatial surface terms read 



dSi $ g v dfi = . (7) 

dG 

There are several ways to achieve the equation above. For Dirichlet boundary conditions one 
would demand $ = at dG and for Neumann type dSig^djQ = at dG. But there is 
also a third possibility for the disappearance of the surface terms, namely if the components 
of the metric themselves which are orthogonal to the surface dSi approach zero at dG, i.e. 
dSi g^ (dG) = 0. As stated above this might be the case for a horizon situated at the boundary 
dG. 
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Strictly speaking, there exist various definitions of a horizon, such as the event, apparent, 



Cauchy, particle, and putative horizon, cf. [33] and [27|. The definition of the particle horizon 



refers to a special observer at a given time-like word-line whereas the other horizons can be 
defined in an observer-independent way. Hence, the vanishing of the spatial surface terms 
(without constraints on the fields) implies a priori only a particle horizon at dG. However, 
with some additional requirements on the space-time, for instance spherical symmetry and 
asymptotic flatness, the (particle) horizon at dG meets the other definitions as well. 
After having established the properties of the space-time M, we arrive at the conclusion that 
it is globally hyperbolic (i.e. strongly causal and complete, cf. [Q) and the spatial domain G 
represents for every time t a Cauchy surface. 

2.2 Energy 

For a time-independent metric the Noether theorem demands the existence of a conserved 
energy. The energy-momentum tensor for the scalar field reads 

t„ = ^= = (0 M *)(a,*) - g -f (9 7 $)(^$) . (8) 

By virtue of the Klein-Fock- Gordon equation □<£> = the covariant divergence of the energy- 
momentum tensor vanishes 

= d, (V=9 TV) - \ T^d u g aP = 0. (9) 

In general this covariant equation does not lead to any conserved quantities due to the exchange 
of energy and momentum between the gravitational and the scalar field (second term). But for 
a stationary metric (<9 g a /3 = 0) it is possible to construct a conserved energy flux utilising 
the [y = 0)-components 

drj" = d ll (y/=jT£)=0. (10) 

This local conservation law allows for the introduction of a conserved energy as a global quantity 
via 

E d ^ J A^T ° = J d 3 rT°. (11) 

G G 

For a Minkowski space-time where T 00 = T° = T 0Q holds this definition coincides (of course) 
with the usual energy. Another argument for the above defined energy for being the correct 
choice is the following: Starting from the action A we may define the Lagrange function L such 
that 

A= [ dtL (12) 
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holds. The Hamilton function H as the Legendre transform of this Lagrange function exactly 
coincides with the energy of the field H = E. 



3 Particle definition 



To provide a canonical definition of particles one has to indicate which properties the particles 
should exhibit. For a free (linear) field we expect the particles to evolve independently and 
to carry a certain energy. As shown in Section |2.2| , for a static metric the energy E of the 
field $ and its Hamilton function H coincide. Consequently, both requirements can be satisfied 
by the diagonalisation of E = H or, equivalently, the Lagrange function L. Having defined 
the particles via diagonalisation of H = E, the corresponding vacuum |0) coincides with the 
ground state of the Hamiltonian and the energy. Of course, the procedure described above does 
not represent the only one possibility to accomplish the particle definition. Another approach 
is based on the "one-particle structure" of classical solutions of the field equation, see e.g. 
H |34|, |7|, [14|, and the remarks in Sec. [3J; . 



According to the definition in Section p.2| the Lagrange function governing the dynamics of the 
field reads 

L= l -j £rg™(r)&+ 1 - j d 3 r ^(r)^)^) , (13) 

G G 

where d 3 r denotes the spatial integration with the volume element d 3 r = \/—g d 3 x. 

To diagonalise this expression one has to deal with an elliptic partial differential operator which 



requires some functional analysis. All of the used theorems can be found in |35[| and are 
not cited explicitly in the following. 



3.1 Hilbert space theory 

To work with mathematically well-defined quantities we have to set up some definitions. C^°(G) 
denotes the set of all infinitely differentiable functions u : G —>■ K of compact support inside 
the open domain G. For two functions of this kind u,v G Cq°(G) we define a scalar product 
via 

{u\v} 1 = J d?r g 00 (r) u*(r) v(r) . (14) 

G 

The assumption of a non-degenerated signature in Section [2.1| is essential for this definition. 
Without a positive g 00 the above expression would be a pseudo-scalar product instead of a scalar 
product with = <-> u = 0. The latter property is necessary for investigations concerning 
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convergence. As every scalar product induces a norm ||m|| 2 = it is now possible to define 

a Hilbert space as the completion of all C£°(G) functions with respect to this norm 

L 2 (G,g°°)^C^{G) { ' l ' } \ (15) 

Because every C ,o (G)-function can be L 2 (G, g 00 )-approximated by linear combinations of step 
functions, this Hilbert space is separable. 

The same procedure may be performed for vector- valued functions u : G —>■ M 3 . Again we 
may define a scalar product for two smooth functions of compact support u, v G [C^^G)} 3 due 
to the non-degenerated signature 

{u\v} 3 d ^ f - f d 3 r^'(r)<(r)^(r) , (16) 



G 



and in analogy the corresponding Hilbert space reads 

LKG,^)^ [C^(G)f { ' l ' h . (17) 



The advantage of the scalar products defined in such a way becomes evident if we use the linear 
partial differential operator 

V: C™(G)cL 2 (G,g 00 ) - L\(G,g ij ) 

<f>(r) - (c^(r)). (18) 

to cast the Lagrange function into the simple form 



1 = \{*\*} 1 -\{ V *\ V *}*- (19) 



Nevertheless, this is still not a representation which is suitable for diagonalisation. For that 
purpose we have to perform the spatial integration by parts (see Section |2|). In terms of 
functional analysis this means the construction of the adjoint operator. The domain of definition 
Def(P) = C^°(G) of the P-operator is dense in L2(G,g 00 ). As a consequence, its adjoint T>' 
exists as a linear operator Z>t : Def(D t ) C Ll(G,g ij ) -> L 2 (G,g 00 ). For [C 3O (G)] 3 -functions 
the spatial integration by parts is always possible. Accordingly, the domain of definition of 
the adjoint T>^ contains these functions [C X) (G')] 3 C Def(P^) an d is thereby also dense in 
L\(G,g l i). Therefore the twice adjoint P'T exists as a linear operator as well P'T : Def(P'^ q 
L 2 (G, g 00 ) — > L\{G, g^). Of course, these operators describe physical reality only if one ensures 



the possibility of the spatial integration by parts via physical reasons as done in Section |2Tl . 
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3.2 /C-operator 

Now we are in the position to cast the Lagrange function into a form which can be utilised for 
the diagonalisation of the system. With the definition of the elliptic partial differential operator 
(see also |TJ, we arrive at 



= liM^^lm^mi- (20) 

Every linear operator of the form K = T>^T>^ is non-negative and self-adjoint and thus can be 
diagonalised. Let us study the domain of definition Def(/C) of this operator. The twice adjoint 
operator is the closure of the original operator D, i.e. D = T>" . Its domain of definition is 
the completion of all C^ (G)-functions 

■{■I-} 



Def (V) =C%°(G) 15 (21) 

with respect to the graph scalar product which is defined via 

{u\v} v = {u\v} 1 + {Vu\Vv} 3 . (22) 

One observes that the operator 1 + /C is exactly the Friedrich extension (which is self-adjoint, 
see fl35l ) of the original operator 1 + /C \cg°(G) mediated via the graph scalar product. As a 
result, if the domain G has boundaries dG with Dirichlet boundary conditions, these boundary 
conditions are already incorporated into the domain of definition of the operators T> and /C, i.e. 

(f) G Def (/C) C Def (V) <f>(dG) = . (23) 

To incorporate Neumann boundary conditions one has to start with an operator like (u 1 ) 1 — > diU 1 
and to proceed in the same way. 

As mentioned in Sec. |^, additional potential terms do not alter the main conclusions. If we 
assume the scalar curvature to be a bounded — m 2 < 71 < and smooth 71 G C°°(G) 

function we may introduce a new operator via 

B: L 2 (G,g 00 ) ^ L 2 (G,g 00 ) 

<f)(r) -> (m 2 +K)<j)(r). (24) 

Obviously this operator is bounded, non-negative, and self-adjoint. In addition, since 
DefiV^V) C Def(£>) = L 2 (G,g 00 ), we may define a modified /C-operator via 

)C^VW + B\ De{(vm , (25) 

which is still self-adjoint and non-negative. 
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3.3 Spectral theory 

As mentioned above, every self-adjoint operator can be diagonalised. One way to reveal this 
statement in a more explicit form is the following theorem: For every self-adjoint operator K, 
there exists a spectral family S of orthogonal projections with 

K = J A d£{\) . (26) 

d£(X) contributes only for values A being in the spectrum a(JC) of the /C-operator A G cr(/C). The 
spectrum cr(/C) of an operator JC contains all complex numbers z G C for which the resolvent 
TZ(z) = f (z — /C) _1 does not exist, i.e. (z — K,)~ l is not a well and densely defined and bounded 
operator. For a self-adjoint and non-negative operator K the spectrum is purely real and non- 
negative cr(/C) C R + . It splits up into two parts, the point spectrum a p and the continuous 
spectrum a c . The point spectrum is the set of all proper eigenvalues A corresponding to proper 
eigenfunctions 

a p = {AGC : 3\f x } : fC \f\} = \f\} A} . (27) 

The continuous spectrum contains all numbers A where (A — /C) _1 formally exists, but is not 
bounded 



{AgCV p : 1 1 (A — /C) 1 1 1 =00} . (2* 



The discrete spectrum o"d is that part of the point spectrum a p which incorporates all iso- 
lated points A of <7 p with a finite number of corresponding eigenfunctions |/a}- The continuous 
spectrum a c may also be divided into two parts, the absolute continuous spectrum <j ac , where 
dS(X)/d\ exists as a weakly integrable operator, and the remaining singular continuous spec- 
trum a sc . 

To provide some physical insight into these abstract quantities we shall investigate the spectrum 
for a few examples. The discrete spectrum <7d describes localised states, such as bound states 
or states of a field confined in a finite volume. The point spectrum a p may contain more points 
with additional characteristics. E.g., if the operator K, governs the dynamics of the Maxwell field 
there is an infinite set of eigenfunctions at the point A = 0. These functions correspond 
to the gauge invariance of this theory and do not change physical quantities. The absolute 
continuous spectrum cr ac represents usually the scattering states, but the singular continuous 
spectrum <t sc may be related to more strange phenomena, like quasi-bound states, scattering 
states in average, fractal measure dfi p (cf. Sec. [J]3] below), chaotic behaviour, etc. 
Fortunately, for smooth and regular coefficients g^ u with an appropriate asymptotic behaviour 
the spectrum of the /C-operator is either purely discrete cr(/C) = o"d (for a finite volume) or 
absolute continuous cr(K) = cr ac (for an infinite volume, see e.g. p4| ). 
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3.4 Spectral theorem 



For our main intention, the diagonalisation of the Lagrange function, it is suitable to make use 
of the following theorem: For every self-adjoint operator K, acting on a separable Hilbert space 
there exists a unitary transformation IA which diagonalises it : UfCW — M.. M. denotes the 
multiplication by argument: (Aif)(\) = A/(A). Because fC is C-real, i.e. (/C0)* = /C(0*), we 
may construct a quasi-unitary transformation 

v= ) : L2(G '"° 0) ~" ©^HK),/v) = ^v), (29) 

which is C-real (V<f>)* = V(0*) and does also diagonalise the operator V/CV^ = M.. 
Accordingly, the Hilbert space L 2 (cr, V) is restricted to real numbers and the associated scalar 
product reads 



{u\v} a = {VMV^K = E / d ^ X ) MA) w p (A) = ^ urw . (30) 

p <r(C) r 

Because ^2(0", V) is a real Hilbert space over K., the usual complex conjugation of the first 
argument in the scalar product disappears. 

For a discrete spectrum cr(/C) = a& the measure g?/i p (A) denotes simply a sum and for an 
absolute continuous spectrum cr(/C) = cr ac an elementary integral possibly together with a p- 
summation, cf. Eq. (pO]). For example, the jp-sum may describe the angular quantum numbers 
for the Laplacian in spherical coordinates p = £,m. Both, the p-summation and the integration 
with the measure d[i p (X) are now abbreviated by the index V. 

Performing the transformation of the fields \Q} a = V |$} 1 the Lagrange function can be diag- 
onalised 



L = \ {Q\Q} a ~\{Q\M \Q} a = ~ ^ ((ft - ulQl) , (31) 

r 

with cjp = f Ar G ^(JC) C ]R + which will be called eigenfrequencies. 

One should note that the \Q} a still depend on time \Q{t)} a1 only the spatial dependence is 
transformed by V. Owing to the reality of the transformation V the amplitudes Qr{t) are real 
as well. 
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3.5 Canonical quantisation 

Starting with the diagonal Lagrange function in Eq. (|31|) we are able to perform the canonical 
quantisation procedure by imposing the usual equal time commutation relations 



[{u\Q(t)}^{P(t)\v 
u\Q(t)\ ,\Q(t)\v 



i{u\ v }„ , 
u\P(t) 



P(t)\v 







(32) 



which hold for all \u] a and \v} a . In this representation the canonical conjugated momenta are 
simply determined by \P} a = \dQ/dt} a . 

Due to the isometry of the transformation V these commutation relations are completely equiv- 
alent to the corresponding relations for the field $. For a static metric the inner product is 
related to the scalar product via 



(V#) = i 



(33) 

As a consequence, the relations above are indeed identical to the commutators of the fields 

v#) , (%)] = m) ■ (34) 

The Hamiltonian splits up into an infinite set of commuting parts describing harmonic oscillators 
that are appropriate for a particle definition 



H=- t P\P 



M 



Q 



(35) 



In terms of the creators 



A* 1 } and annihilators 



i} = (M 1/4 = 0)} + iM~ 1/A \P(t = 0)} 



(36) 



the Hamiltonian can be cast into the form 



H 



Ai 



M 1/2 



A 



1 

2 



{A\M^\ti} r ^ (4i r + i r 4 



(37) 



For a discrete spectrum a{K) = ad this already defines the physical particles because we have 
now creation and annihilation operators A\, and A? that diagonalise the Hamiltonian (which 
is also the energy operator). 
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For a continuous spectrum er(/C) = a c the quantities Pp, Qp, Ap and A r are not well-defined 
operators but operator-valued distributions because 

Q T (t),P A {tj\ =i6(T,A), (38) 

which is a Dirac 5-distribution for continuous indices T, A. But a product of two distributions 
acting on the same linear space (e.g. the Schwartz/Sobolev space &i for the 5-distribution), 
i.e. with the same index T, is not well-defined. This reflects the infinite-volume divergence in 
quantum field theory. Consequently, the Hamiltonian in Eqs. (|35|) and (|3~?D is not well-defined. 
It may only be viewed as a formal expression until an appropriate regularisation method has 
been applied. 

3.6 Vacuum definition 

In order to get rid of the singularities discussed above and to obtain well-defined operators dj 
we introduce a complete orthonormal and real basis \bj} a with I G N of the separable Hilbert 
space L 2 (cr, V) and define 

a/^jiifej} . (39) 

For a discrete spectrum cr(/C) = era we may choose bj(T) = Spi which leads us back to the 
operators Ap. For a continuous spectrum er(/C) = a c this coincidence does not hold. Due to 
{bi\bj} a = w ith a Kronecker-5/j the di are well-defined operators with [di, dj] = Su in- 
stead of operator- valued distributions with [Ap, A\] = S(T,A). Unfortunately, for a continuous 
spectrum cr(K) = cr c the operators di are now well-defined, but do not exactly diagonalise the 



Hamiltonian. But - as we shall see later in Section |3.7| - one may choose an appropriate basis 
\bi} a for which the operators di approximately diagonalise the Hamiltonian. 
The corresponding number operators take the usual form hi == d\dj. The Fock space $ which 
contains all pure states of the quantum field $ is now defined as the completion of the 
linear hull of the proper eigenvectors of these commuting operators hi for all indices I G N 



5 = lin ) : V/ : n/ |tf ) = n/} . (40) 

As a consequence, the spectrum of the operators hi in this Fock space is a pure point spectrum 
o~(hi) = <t p . With the same arguments as already used for the quantisation of the harmonic 
oscillator the commutation relations [a/,Oj] = Su imply er(nj) = N. It should be mentioned 
here that this definition of the Fock space $ is slightly different to the frequently employed 
approach based on the one-particle Hilbert space Sj ~ L2 (G, g 00 ) (see e.g. |], 0, 0, |l^, |15| ) 

J = £(BS) © (i3 ©fl) symm © (3 ©^ ©-^) symm © • • • • (41) 
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Nevertheless, these distinct definitions are related if one divides the Fock space in Eq. (f40|) 
into orthogonal subspaces labelled by different values of the total number of particles n total = 

Accordingly, the vacuum is defined as the eigenvector of all commuting operators "hi with 
eigenvalue zero 

V/ : hi\0) = i.e. ctj |0) = 0. (42) 

This definition is independent of the special choice of the basis \bi} a . To prove this statement, 
we use the completeness of the basis \bi} a to obtain the following result 



V ICL : (A|CjJO> = 0. (43) 

If we regularise the formal expression for the Hamiltonian in Eq. ([3^]) via insertion of a complete 
basis (in principle together with a convergence factor, exp(— Ie) for instance) 

H = {{M^A^h] {b!\M 1/4 A} +{.M 1 / 4 i|& / } [bjlM^A^ ) 

= Y^^\M^%} {m^%\A) + l -{M l /%\M l '%}^j , (44) 

it appears as a divergent sum of some non-negative operators of the structure XjJtj and re- 
maining C-numbers. The above defined vacuum is the ground state of all operators XjXj 
and in this regard also the ground state of the Hamiltonian. Hence, the divergent amount of 
C-number terms represents the zero-point energy. The infinite summation over the index / 
corresponds to the sum over arbitrary high frequencies and - for a continuous spectrum - the 
summation of an infinite number of basis elements for a given frequency interval. The first 
infinity, the infinite energy divergence, is always present in quantum field theory and the latter, 
the infinite volume divergence, only for non-discrete spectra. 

In the Minkowski space-time the above defined vacuum coincides (of course) with the usual 
Minkowski vacuum |0) = |0m)- In the Schwarzschild space-time this state - which is the ground 
state of the Hamiltonian - is called the Boulware f| state |0) = |^b)- 

The particle definition presented above can be reproduced utilising the well-known approach 
based on the inner product: The basis elements \bi} a of the Hilbert space L 2 (cr, V) are nor- 
malised and therefore correspond to functions ej(r) via |e/} 1 = \bi} a which are also nor- 
malised {e/lej}^ = 6u and build up a basis of the Hilbert space L 2 (G,g 00 ). As a consequence, 
the operators aj correspond to localised wave-packets \F I (t)} 1 which are defined as follows 

= (4/C)" 1/4 exp {-iK}'H) Ie,}! . (45) 
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These quantities are solutions of the Klein-Fock-Gordon equation dj \F I } 1 = —JC \F I } 1 and 
normalised with respect to the inner product 

(F i \Fj) = -(F*\F})=5 I j , (F*\Fj) = (Fj\F})=0. (46) 

Comparison with the particle definition via the inner product verifies indeed the identification 

a 7 = . (47) 

The functions Fj and Ff form a complete set of solutions of the Klein-Fock-Gordon equation. 
Hence, the field $ may be expanded via 

$ = ^^ + 4^*, (48) 
I 

which demonstrates again the equivalence of the approaches. 
3.7 Eigenfunctions 

For a point spectrum a p there exist proper eigenfunctions fr G L 2 (G,g 00 ) with K,\fr} 1 = 
u>y \fr}i, but for a continuous spectrum a c this is of course not the case. Nevertheless, it is 
in many cases possible to find an analogue. If pointwise defined functions fr(r) (or - more 
generally - locally integrable functions fr G L l ° ca - 1 ) exist such that 

{C|V0} CT = ^Cr| rf 3 r/ r (r)0(r) (49) 
r G 

holds for all <fi G C^°(G) and ( G C£°(cr) the functions /r( r> ) are called (generalised) eigen- 
functions of the /C-operator. In contrast to the proper (cr p ) eigenfunctions with fr G Def(/C) C 
L 2 (G, (7 00 ) the generalised (<r c ) eigenfunctions do not belong to the Hilbert space fr G" L 2 (G, g 00 ) 
and (of course) also not to the domain of definition of the /C-operator / r G" Def(/C). However, 
due to VIC — MV also the generalised eigenfunctions fulfil the pointwise/local (generalised) 
eigenvalue equation /Ci oca i/r(^) = ^r/r(^)- This is a very important relation for the calculation 
of these eigenfunctions. If the (generalised) eigenfunctions exist, the transformation of the fields 
|$} 1 = V" 1 " \Q} a can be described by the pointwise/local identity 

Hr,t)= ^QrWMr). (50) 
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Even though the generalised eigenfunctions are not in L 2 (G,g m ) ) they may be thought as 
a (singular) limiting case of L 2 (G, g 00 )-functions: In the following considerations we assume 
dfj,(\) = d\ and a = IR for reasons of simplicity. The L 2 (cr, V)-basis functions bj(\) can 
be squeezed and translated b £ j^(\) = b I (\/e — \r/e)/y/e and are still a basis of L 2 ((J, V). 
Evaluating the (singular) limiting case of these squeezed basis functions 

lim = lim MA^-ArO = _ 

where A//(r) denotes some normalisation factor, we observe that every generalised eigenfunction 
fr(f) can be locally approximated by appropriately chosen wave packets ef(r)| = ^ ^f(r)} 



Iim-ffi4=Mcr)/r(r). (52) 

£j.u y £ 



Accordingly, also the operator- valued distributions A r may be considered as a singular limiting 
case of the regular operators 



a 6 , 



lim = A/} (r) i r . (53) 

The divergent factor 1/y/e indicates the singular character of the generalised eigenfunctions 
(e.g. plane waves) in contrast to the regular basis elements (wave packets). Of course, in 
realistic experiments one never deals with plane waves, but wave packets. On the other hand, 
the calculations with plane waves are usually much simpler. Hence, in the following we shall 
perform our evaluations with eigenfunctions always bearing in mind their character as a singular 
limiting case of regular objects. 



3.8 Continuum normalisation 

To investigate the physical consequences caused by the singular behaviour of the product of two 
distributions iVr = i r ir - expressed by the factor \ je - we consider a quantum field confined 
in a finite volume V and study the limiting case V — > oo. This limit may be interpreted as 
the transition from a discrete spectrum cr(/C) = ffj to a continuous one a(/C) = a c . For a 
3-dimensional cubic volume V the indices T correspond, for example, to discrete wave-numbers 
k. In the continuum limit V — > oo the fc-sum transforms into an integral over d 3 k via 




Afy denotes a normalisation factor which depends on the imposed boundary conditions (Dirich- 
let, Neumann, periodic, etc.) and the shape of the domain G. The Kronecker-<5i, y converts 
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into a Dirac-<5 3 (A; — k') in an analogue way NvVSuu' —> $ 3 (k — k'). Ergo, the singularity 
5 3 (k — k) displays the infinite-volume divergence 5(k,k) = 5 3 (k — k) = MyV . Recalling the 
formal expression for the Hamiltonian in Eq. (|57|) 



it J , ! f /p ^ ^ A. ,jv , 



if = ^r( AMr + 2 5 ( r ' r ) J = £ + -±6{T, T)j , (55) 

r 

we observe that - in addition to the mode summation/integration - its indefinite character 
exactly exhibits this divergence. Indeed, if one examines the continuum limit of the Hamiltonian 

kA = E (H + 1) — > ^ c = / rf3fc (^fe + ^ 3 ( fc - ' ( 56 ) 

the limiting number "operator" can be identified viaAfyVN^ — > Ni. The singular charac- 
ter of this formal expression may be exemplified with the following consideration: If the state of 
the quantum field corresponds to thermal equilibrium at a temperature T > 0, the expectation 
value of the number operator equals the Bose-Einstein distribution for arbitrary large but 
finite volumes V. For an infinite volume the expectation value of the quantity diverges 
owing to the factor MyV. The expectation values of the regular operators ni are (of course) 
still finite and behave as the Bose-Einstein distribution evaluated at some averaged frequency 

Ui. 

3.9 Bogoliubov coefficients 

So far we have considered static space-times and developed an appropriate particle definition. If 
we now drop the restriction to stationary metrics and take dynamical space-times into account, 
the question concerning particle creation arises. A variation of the metric g^ v induces a change of 
the /C-operator and - possibly - the corresponding Hilbert space L 2 (G, g 00 ). A function, which 
belongs initially to L2(G m , gf®) may be later (e.g. if a horizon has formed) not in L 2 (G out , g^ut) 
but a distribution with respect to the L 2 (G out , <j£° t )-scalar product. As a consequence, it is not 



clear whether the Bogoliubov coefficient (see e.g. p^]-|5l|) describing the particle creation 

P u = (F;\Fj) (57) 

exists for all Fj G L 2 (G m , gf®) and Fj G L 2 (G out , g°® t ) or not. However, C^-functions belong 
to the domain of definition of all (tempered) distributions. Thus for Ff 1 G C^°(G m ) and 
Fj ut G C^G 011 *) the Bogoliubov coefficients always exist provided the metric can be cast 
into an analytic form. Similar to the previous Sections all other quantities (e.g. generalised 
eigensolutions) have to be approximated with C^-functions. 
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4 Flat space-time examples 



In the previous Section we have derived a canonical definition of particles for curved space- 
times which fulfil certain conditions. In the following we are going to apply this approach 
to the most simple example of a flat space-time in order to achieve a deeper insight into the 
physical consequences of the used mathematical theorems. 

For the unbounded 1+1 dimensional Minkowski space-time with ds 2 = dt 2 — dx 2 the /C-operator 
reads 

(»> 

together with the domain G = (— oo < x < oo). The infinite volume of this domain and the 
regularity of the metric cause a purely absolute continuous spectrum cr(/C) = cr ac = R+. The 
unitary transformation U (see Section |3.4| ) is simply the one-dimensional Fourier transformation 
U = T . The quasi-unitary transformation V takes the real and imaginary parts separately 
leading to the generalised eigenf unctions sin tux and cos tux. Hence, the expansion of the field 
l> takes the following form 



t) = ^ Qu,c{t) cos(tux) + Q W p(t) sin(tux) 



/7 
77= (AL, c e*"* cos(tux) + Al ;S e iwt sin(tux) + h.c.J , (59) 

o 

with a normalisation factor A/^ depending on the explicit form of the measure d/i(tu 2 ), e.g. 
dfi(u 2 ) = du or dfi{uj 2 ) = dui /2n etc. The spectrum of the /C-operator discussed above is twice 
degenerated, i.e. there are two independent generalised eigenfunctions (sin tux and cos tux) for 
every point A = tu 2 of the spectrum a . This degeneracy of the spectrum allows for the definition 
of particles with a definite direction of propagation: With a simple linear transformation we 
may rearrange the expansion of the field 

oo 

/A 
(Al t+ e iut+iu)X + Al t _ e iwt ~ iwx + h.c.J . (60) 

o 

The new introduced quantities A Wt+ , AL+ and Al > _ obey the same commutation relations 

as the original ones A W)C , A ws , and A* s . Thus they also describe particles. In contrast 
to the original particles which correspond to standing waves with different phases (sin tux and 
cos tux) the new particles describe left-moving and right-moving waves according to exp(±ztux). 
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These complex functions are suitable for a definition of particles A? but do not correspond to 
Hermitian amplitudes Qr- 

As a second example we study the situation of a bounded domain G = (0 < x < oo) in a 1+1 
dimensional Minkowski space-time with a Dirichlet boundary condition (a mirror) at x = 0. 
Even though JC seems to have the same form as in Eq. ( |58l) it denotes a different operator as 
a result of the boundary condition. With the same arguments the spectrum is purely absolute 
continuous o = cr ac = R + . In contrast to the previous example this spectrum is not degenerated. 
Every point A = uo 2 of o corresponds to exactly one generalised eigenfunction, i.e. sin tux. As 
a consequence, the definition of particles with a certain direction (left-moving or right-moving) 
is not possible. This result is physical reasonable if one takes conservation laws into account. 
Every left-moving component will be reflected by the mirror at x = after some period of time 
and turns its direction into right-moving and vice versa. 

A finite domain G=(0<x<L)inal+l dimensional Minkowski space-time with Dirichlet 
boundary conditions at x = and x = L of course possesses a purely discrete spectrum 
a = 0d with proper eigenf unctions ~ sin(N7ra:/L). The insertion of mirrors represented by 
Dirichlet boundary conditions usually lowers the "density" of the spectrum a, i.e. the number 
of eigenfunctions. 



4.1 Ingoing and outgoing particles 

Now we shall extend our investigations to the 3+1 dimensional Minkowski space-time described 
by different coordinate systems. Using spherical coordinates r, <p it will turn out that the 
definition of ingoing or outgoing particles is not possible within the canonical approach. This 
is a consequence of the spectral properties of the operator 

K = -V 2 = -|^, (61) 

together with the domain G = M. 3 . Expressed by Cartesian coordinates r = (x, y, z) T the 
generalised eigenfunctions take the simple form sin(fcr) and cos(fcr) with \k\ = u. As it is 
well-known these functions form a complete basis of /^(K 3 ). 

Employing spherical coordinates r, (p the Cartesian eigenfunctions can be expanded with the 
aid of the equality 

oo 

exp(ifcr) = i£ ( 2i + l)je{ur)P £ {cos 9) , (62) 
e=o 

with kr = lot cos 6 and the Legendre polynomials P%. By inspection we recognise the following 
fact: For a given angular behaviour the spherical Bessel functions je(ujr) are already complete 
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to describe the radial dependence. The Neumann n^(wr) or Hankel functions hf(ur) are not 
required and would be "over-complete". Therefore they do not describe additional degrees of 
freedom and do not enter the particle definition. This result can also be derived by considering 
the spectrum of the /C-operator. Due to the singular behaviour of the functions n^ur) and 
hf(u;r) at r = they are not eigenfunctions. 

To acquire real eigenfunctions we have to introduce redefined spherical harmonics ^m.(^, f) == 
A^ m (costf)cos(ra</?) for m > and ^ m (tf , </?) = A^ m (costf) sin(m^) for m < 0. de- 
note the associated Legendre polynomials and Afe are normalisation factors. Accordingly, the 
complete set of real and orthogonal eigenfunctions reads 

Ue m (r, <p) = N u t3t{wr)yt m {&, y?) . (63) 

Again we observe the occurrence of exactly one eigenfunction f^im per eigenfrequency uj for a 
fixed angular dependence £, m. The regularity at r = plays the role of an effective boundary 
condition and forbids the existence of additional eigenfunctions such as n^ujr) or hf(ujr). As a 
consequence, within the canonical approach it is not possible to define radial ingoing or outgoing 
particles in the Minkowski space-time. Functions like exp(±iujr)/r are not eigenfunctions of 
the Laplacian and therefore not solutions of the wave equation 

B eX ^ Wt r ±tUJr) =-Ane^(r)^0. (64) 

Expanding the field $ into functions that do not satisfy the equation of motion □$ = would 
abandon the independence of the distinct particles. Functions like exp(±iujr) /r correspond to 
the resolvents TZ(uj 2 ± is) of the operator /C (remember <r(/C) C M). Particles are defined with 
respect to the eigenfunctions which are representations of the spectral family £ of the operator 
/C. Into this spectral family S the resolvents themselves do not enter, but linear combinations 
of them: 7Z(\ — is) — 1Z{\ + is) — > £(A) which again leads to sin(cjr)/r = uj ji =Q (ur). 
The impossibility of defining radial ingoing and outgoing particles is not restricted to the 
Minkowski space-time, this holds also for arbitrary spherically symmetric metrics 

ds 2 = g 00 (r) dt 2 + g u (r) dr 2 + r 2 dtt 2 , (65) 

provided that the coefficients of the metric g 00 and g u are smooth and analytic functions. Such 
functions can be Taylor expanded 

2 

9oo(r) = <?oo(0) + <4(0) T - + G(r 3 ) , (66) 

where g' 0Q (0) and ^^(O) have to vanish for smoothness. After the separation of the angular 
variables with ye m ($, f) the radial dependence of the eigenfunctions is governed by a second- 
order ordinary differential equation in r. Provided its coefficients are smooth and regular 
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the solutions of such an equation are uniquely determined by the first two (non-vanishing) 
terms of their Laurent expansion. For the evaluation of these initial data only the terms 
#oo(0), #00(0) — 0, (711(0) and g' u (0) = are of relevance. For that purpose the radial part 
of the metric can be approximated by ds 2 = goo(0) dt 2 + gn(0)dr 2 . Ergo the behaviour of 
the corresponding eigenf unctions is (up to a simple scale transformation with g 00 (0) and gn(0) 
respectively) asymptotically (r — > 0) the same as in the Minkowski space-time. Consequently, 
also in these more general spherically symmetric metrics there exists exactly one eigenfunction 
for given u, £, m which forbids the definition of radial ingoing and outgoing particles. 
In view of conservation law arguments the nonexistence of ingoing and outgoing particles in 
regular space-times appears very plausible: Every ingoing component will bounce off at the 
origin after some period of time and eventually turn into outgoing and vice versa. 



4.2 Rindler metric 



As stated in Section [2.i| the particle interpretation crucially depends on the selection of a 
particular time-like Killing vector. In the following we shall consider an example where this 
dependence will become more evident. In the previous treatments we focused on the Killing 
vector mediating the Minkowski time translation symmetry. Of course this Killing field cor- 
responds to usual observers at rest. But there exist further time-like Killing vectors in the 
Minkowski space-time - associated with special Lorentz boosts - which result in a deviating 
particle interpretation. 

Starting with the 1+1 dimensional Minkowski metric ds 2 = dt 2 — dx 2 and performing the 
coordinate transformation 

t = p sinh kt , 

x = p cosh kt , (67) 
one arrives at the Rindler metric ds 2 = K 2 p 2 dr 2 — dp 2 . The quantity k is called the surface 



gravity, see e.g. p8fl . For fixed p the transformation describes an accelerated motion. With 
respect to the new time coordinate r the Rindler metric is static and thus allows for a particle 
definition according to Section ^. These particles may be interpreted as those seen by an 
accelerated observer. The corresponding /C-operator can be cast into the form 

*: = -«»„!,!, ( 68) 

with G = (0 < p < 00). The surface term (see Sec. |2|) at p = vanishes without imposing any 
condition on the field $ due to yj—g = up = at p = 0. Indeed, the Rindler metric possesses 
a horizon there. Since the occurrence of this horizon depends on the choice of the coordinates 
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and thereby on the observer, it is a particle horizon (with respect to all world-lines p = const) 
but not an event (or apparent, etc.) horizon, see Sec. ^J] and [^7], j33fl . 

For further investigations it is convenient to introduce the tortoise coordinate p* = ln(/tp)//t. 
In terms of this coordinate the metric reads ds 2 = e 2Kp * (dr 2 — dpi) resulting in the operator 



with G = (—00 < p* < 00). As a consequence, the spectrum is twice degenerated and the cor- 
responding eigenf unctions read sin^p*) and cos(u;p*) or exp(±iup*), respectively. Returning 
to the coordinate p the eigenf unctions behave as exp(±iu\ii(Kp)/ k) = (np) ±tul ^ K . Even though 
the domain is bounded G = (0 < p < 00) there are two eigenf unctions per eigenvalue which 
allows for the definition of left-moving and right-moving particles. This indicates the absence 
of real boundary conditions on the field <3? at the horizon p = 0. In this regard the horizon 
is the opposite of a mirror. Even for a finite domain G = (0 < p < L) the spectrum of the 
/C-operator is still continuous due to the horizon: G = (—00 < p* < L*). 

4.3 Unruh effect 

After having performed a particle definition for the Minkowski and the Rindler observer, the 
question about the relationship of these two approaches arises. Evaluating the expectation value 
of the number of Rindler particles in the Minkowski vacuum one obtains a thermal distribution 
function, a consequence of the Unruh 0] effect. This effect demonstrates manifestly that 
different observers may obey distinct particle interpretations. In consequence the vacuum may 
depend on the particular Killing vector. 

One way (see also []3T|) to calculate the expectation values explicitly is based on the Bogoliubov 



coefficients 

Ptk = i J F r M dp Ff. (70) 
The generalised Minkowski eigenfunctions are labelled by Y = (£, uj) 

F^(x) = t) = AV eXp( ^ t} c*" , (71) 

- 'UJ 



where £ = ±1 distinguishes the left-moving and right-moving particles. {Mm denotes a normal- 
isation factor.) In analogy the generalised Rindler eigenfunctions read 

Ff{x) = if v , )( p, r) = ^M-^'r) M gv /W (72) 
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With the choice for the surface E = {r = 0,0<p< 00} the surface element takes the form 
d~E° = g 00 y/—g dp = dp/ {up). At this surface the Minkowski coordinates are simply given by 
t = 0, x = p and the derivative transforms according to d T = npd t . Putting all this together, 
the /^-coefficient transforms into 

00 

%«);(£V) = ^MR / , e p {Kp) I . (73) 

J Kp vww' 


This integral involves generalised eigenf unctions (corresponding to A-p) and has to be under- 
stood in a distributional sense. For well-defined expressions (such as di) we have to insert a 
convergence factor, for instance (np) e exp(-EKp). The existence of the limit e I confirms 
the possibility of approximating the singular eigenfunctions by regular quantities. After this 



procedure we may make use of the formula E(d 



due~ uw u z ~ L =w- z T(z), (74) 







which holds for 9?(u>) > and ^St(z) > 0, and we arrive at 

/WeV) = ^mr^ J^- T^'lu'/k) / k + e)~* u ' . (75) 

K V OJ 

Calculating the remaining Bogoliubov coefficient cn^yf^^i) one gets nearly the same expression 
but with a positive sign in front of the term i£u/n. Therefore both coefficients merely contribute 
for particles moving in the same "direction" /% )U /);(£',u;') ~ respectively, a^ iU1 yui iU1 i^ ~ S^g. 
Now it is possible to compare both quantities. As said before, the only difference between a 
and (3 is the sign in front of the term oj/k. Dividing the two coefficients all other terms cancel 
and the convergence factor e determines the side of the branch cut of the logarithm in the 
complex plane. Hence we find 

= exp(— 7T0// k) Ot(Z>,wW,w') ■ (76) 

An alternative way to obtain this important result is based on the analytic continuation into 
the complex plane. For that purpose we define slightly modified Bogoliubov coefficients via 



(5 c u ,(u, u') = \/uu?(3^ u y ){e>u/ ) , (77) 

and in analogy the a-coefficient. In view of Eqs. (fT0|)-(|73]) the modified Bogoliubov coefficients 
can be analytically continued to the complex cu-plane where the relation 

a% ii ,(u,u/)=l3l ) z,(-u,u/) (78) 
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holds. Inserting this equality into Eq. ( [75] ) reproduces Eq. (|76|). In order to evaluate the 
absolute value squared of the /^-coefficient we may utilise the identity ]36 



T(z)T(-z) = — (79) 

Z Sin 7TZ 

to obtain the final result 

■ 1 2 _ 8ttA^ r S u > 

\P(t,u);(e,u')\ 77S 77^ 7- t° U J 

11 huj exp(2TTu'/K) — 1 

In view of the remaining ^-integration the number of Rindler particles in the Minkowski vacuum 
diverges. This result can also be re-derived using the well-known unitarity relation 



«Ar«xr - PatPtt = &{K T) , (81) 
r 

where F symbolises the Minkowski index. This equality reflects the completeness of the 
Minkowski solutions. Special care is required concerning the derivation of an analogue expres- 
sion involving the Rindler functions since these solutions are restricted to the Rindler wedge 
and thereby they are not complete in the full Minkowski space-time. 
Inserting Eq. (|76|) and considering the singular coincidence A = T it follows 



= (0m| ^',w') I°m) = j2 |/%u/);(£',u>')| 2 

6(w,u) AT V V 



exp(27ru//«) — 1 exp(2iru' / k) — 1 



(82) 



In the last step we have used the results of Section fi.8\ Recalling the argumentation made 
there we come to the conclusion that the divergence of N-p is necessary for a thermal behaviour. 
The same calculation can be performed with well-defined operators nj corresponding to lo- 
calised wave packets. For an appropriately chosen basis ej(r) the coefficients (3u are up to 
normalisation factors approximately the same as the (5yk evaluated above. But in this case the 
results for (hj) = J2j \Pu\ 2 are finite owing to 5n = 1. Another explanation is the fact, that 
the flu for arbitrary frequencies u do not coincide with the firA (due to the localised character 
of the wave packets) which makes the (^/-summation finite. One way to perform technically 
such a calculation involving localised quantities is to insert a convergence factor with a finite e 
similar to the comment after Eq. (|73D . Accordingly this finite e enters the /3-coefficients and 
causes a finite result of the tuj-summation and thereby a finite number of created particles as 
well. Omitting the corresponding normalisation factor the infinite volume divergence can be 
restored in the limit e 1 0. 
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4.4 KMS condition 



From a strictly axiomatic point of view the divergent result in Eq. (p2|) in the last section 
may not be completely convincing. However, it is possible to show more rigorously that the 
Minkowski vacuum indeed behaves as a thermal state when analysed by a Rindler observer. 



This can be done by employing the Kubo-Martin-Schwinger (KMS) condition 37, 38]. A 
KMS state (-)t is defined as a time-translationally invariant state which satisfies the following 
condition 

(U(t)V(t')) T = (V(t')U(t + i/T)) T (83) 

for all observables U and V and some temperature T. It can be shown that if the (irreducible) 
algebra of observables possesses a well-defined matrix-representation then the KMS state cor- 
responds to the usual canonical ensemble 

0)T= Ju?M^m\ ( 84 ) 

One might wonder at the fact that the Minkowski vacuum, i.e. a pure state, displays thermal 
features - usually connected with mixed states. This can be explained by the thermo-field 
formalism, see e.g. [B9||. I.e., a pure state of a quantum system transforms into a mixed state 



after averaging over a subsystem owing to the correlations between the different subsystems. 
As a result of the particle horizon at p = the Rindler observer is causally separated from a 
part of the Minkowski space-time and does therefore indeed regard the Minkowski vacuum as 
a mixed state. 

To show that the Minkowski vacuum displays the temperature T = k/(2tt) we consider the 
corresponding two-point Wightman |4TJ function. The Wightman axioms (in particular the 



spectral condition) imply that this bi-distribution can be considered as the boundary value of 
an analytic function. Hence we may restrict to the space-like region for reasons of simplicity 
where the two-point function assumes the form 

W(x,x!) = (*(x)*(x!)) = - / 2?r) ' (85) 

for 3+1 dimensions. In 1+1 dimensions it behaves as ln[(a; — x') 2 ) which does not alter the 
following considerations. Since for a free field all n-point functions can be derived from this 
2-point function it contains all information about the theory. 

Now we may consider the two-point function in terms of Rindler coordinates. The t, x- 
contribution to the geodesic distance transforms according to Eq. (|67|) into 

{t-t' f -(x-x'f = 2pp' cosh (k[t - r ]) - p 2 - p' 2 . (86) 
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As a result the two-point function is periodic along the imaginary Rindler time axis and thus 
satisfies the KMS condition for the temperature T = k/(2tt). This result confirms the consid- 
erations in the previous section and justifies the identification of the UV-divergence occurring 
there with the infinite volume divergence of the Rindler space. 

It can be shown quite generally that only the KMS state corresponding to the temperature 
T = k/(2tt) satisfies the Hadamard [^0| condition (local stability) in the complete Rindler 
space-time (and in particular at the horizon), see [fL3|l . The Hadamard requirement demands 
the singularity of the two-point function (1/s 2 and Ins 2 ) to be independent of the state, i.e. it 
is only determined by the structure of the space-time, see also |17|, |23|, [25|, This property 
ensures the validity of the point-splitting renormalisation technique, cf. |IIJ. As it will become 
more evident in Section BTTL an analogue idea can be employed to derive the Hawking effect. 



5 Black holes 

In this Section we are going to apply the formalism presented in Sec. |3| to one of the most 
fascinating curved space-time structures, the black hole. Various coordinate systems which 
represent this object are known. For our purpose we have to demand a static metric with a 
time coordinate t corresponding to a Killing vector. Because the black hole space-time becomes 
asymptotically flat, another requirement is the coincidence of this time coordinate t with the 
usual Minkowski time of an observer at spatial infinity. All this requisites are fulfilled by the 
Schwarzschild coordinates t, r, ip for which the black hole metric reads 

dv^ 

ds 2 = h(r) dt 2 - - r 2 d® 2 - r 2 sin 2 d^ 2 . (87) 
h(r) 

Other coordinates, e.g. Kruskal, Eddington-Finkelstein, etc. are not suitable for the above 
reasons. As the Schwarzschild coordinates measure time and length scales with respect to 
an observer at fixed spatial distance to the black hole all results obtained later refer to this 
observer. 

In order to describe a (non-extreme) black hole with a horizon at r = R and a surface gravity 



k > the function h obeys the properties (see e.g. |28 



h(R) = , k = - h\R) (88) 

and also h(r > R) > together with h(r — > oo) = 1. With the aid of this function h it is 
possible to consider the rather general case of a static black hole, for example the Schwarzschild 
metric with h — 1 — R/r. 
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Using these coordinates the canonical conjugate momenta turn out to be ft = dt^/h. In 
terms of these momenta the formal expression for the Hamiltonian density can be cast into the 
following form 

* = + £ (*i) 2 + ± ^ (a^) 2 . (89) 

The fields $(r,£) as well as their momenta fl(r,t) are operator-valued distributions (see also 
Section ^). Consequently, the Hamiltonian density above is not well-defined. In analogy to 
Sec. p.5| it may only be considered as a formal expression until an appropriate regularisation 
method, for instance the point-splitting technique (see e.g. flQfl), has been applied. 
It is possible to split up the Hamiltonian H of the field l> into two parts H = H > + H < that 
account for the interior H < and the exterior H> region of the black hole, respectively 

H> = J d 3 r H 6(r - R) , 

H< = I d 3 r H e(R-r) (90) 



with the Heaviside step function and the volume element d 3 r = y/^g d 3 x = r 2 sin •& dr d$ dip. 
Employing the equal time commutation relations 

Mr, t), <f>(r', t)} = [fl(r, t), fl{r', t)} = ; [$(r, t), fl(r', t)} = i5 3 (r - r') , (91) 

where t denotes the Schwarzschild time and represents a Killing vector, one observes that the 
two parts of the Hamiltonian commute 



. (92) 



In the language of point-splitting, cf. |TIJ and the remarks in Section O, the divergent terms of 



the Hamiltonian density are independent of the state and therefore pure C-numbers which do 
not contribute to the commutator. The remaining (convergent) operator-valued components 
commute because of h(r = R) = 0. The same result can be obtained by means of normal 
ordering or the regularisation described in Eq. (f44j). Due to h(r = R) = the /C-operator 
and the operators projecting onto the interior, respectively, exterior domain commute. Hence 
it is possible to select a basis hi for the inside and outside region separately such that the 
Hamiltonian possesses no mixing terms. 

Accordingly, the separation H = H > + H < represents two independent systems. This fact 
displays one advantage of the Schwarzschild coordinates because there is a horizon at r = R. 
The consistency with the results of Section |2J] can be demonstrated if one considers the spatial 
surface term dSi = d 2 x \J —g rti g^ = dd dip r 2 sin d g rr which indeed vanishes for r = R. As 
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a consequence, it is impossible to transport matter (energy or information) across the horizon, 
nothing can come out or fall into the black hole. Of course, this holds only for a fixed metric, 
i.e. if one neglects the back-reaction. Without this restriction it is possible that the horizon 
increases due to the in- falling matter, and swallows it. It should be emphasised again that all 
of our assertions refer to an observer at a fixed spatial distance to the black hole and therefore 
not necessarily to a free falling one. 



5.1 Black hole exterior 

In the following we restrict our considerations to the domain outside the black hole governed 
by H > . The properties of the interior will be discussed in the next Section. The exterior region 



G = {r > R} fulfils the conditions imposed in Sec. [2.1| which allows for a particle definition. 
As a result, the f/>-part of the Hamiltonian can be diagonalised formally via 



#>= T tu r N™ + , (93) 



where Eoo denotes the divergent zero-point energy. 

In order to isolate the features that are specific for black holes, the most interesting region is 
the neighbourhood of the horizon r « R. To investigate the behaviour in this zone we introduce 
a dimensionless variable \ with 

x = 2K(r-R)^h = X (l + 0(x)) ■ (94) 

Without loosing the static character of the metric it is possible to perform a radial coordinate 
transformation for r > R via 

¥ = ^ + 0(x). (95) 

The new radial r* coordinate is called the the Regge- Wheeler tortoise coordinate. According 
to the above arguments it is sufficient to cover the region outside the horizon by the new 
coordinate. The function h and the original radial variable r have to be considered as functions 
of the introduced coordinate: r = r(r#) and h = h(r*) = /i(r[r*]). The tortoise coordinate has 
the advantage of a very simple form of the /C-operator 

K =4sr 2 ^- w ^=-l? +ow ' (96) 
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together with G = {r* G R}. The generalised eigenfunctions F® u (x) of this operator behave 
(9(x)-approximately as exp(±iur*) and after the separation of the angular variables they can 
be written as follows 



F r BH (x) = Ff* m (t, X , 0, <p) = ACT -fc X^ /(2K) y*n(*, <P) (1 + 0(x)) ■ (97) 



A(T symbolises a normalisation factor which may without any loss of generality chosen to be 
independent of £. These eigenfunctions are rapidly oscillating near the horizon. 
By inspection, we recognise the occurrence of two generalised eigenfunctions for a given fre- 
quency us and fixed angular dependence £,m distinguished by £ = ±1. Thus the definition of 
ingoing and outgoing particles is possible in this case. (It should be mentioned that potential 
scattering effects cause slight deviations from the purely ingoing and outgoing behaviour in Eq. 



(p7|) at r — ► oo. However, this way of definition does not alter the conclusions.) This - perhaps 
surprising - fact can be elucidated in the following way. The horizon separates the space into 
two independent domains (interior and exterior) and prevents the field modes outside from 
being influenced by the effective "boundary condition" at r = 0. In view of the study of the 
/C-operator in terms of the tortoise coordinate r* one may consider the horizon as some new 
kind of spatial infinity (r* — > — oo) in addition to r — > oo. 

Also for the black hole example the horizon acts opposite to a mirror, cf. Section |4.2j . Even for 
the scenario of a black hole which is enclosed in a large box with Dirichlet boundary conditions 
the spectrum of the operator K is still continuous - but now not degenerated. 
For a black hole in an asymptotically flat (unbounded) space-time there are two contributions 
to the infinite volume divergence (see Sec. |3.8|) 5(T, T) = AfyV: firstly, the usual infinity 
r, r* — > oo and secondly, the effective infinity at the horizon r — > R respectively r* — > — oo. The 
former divergence 8 + (T, T) does also arise in the (unbounded) Minkowski space-time - but not 
inside a finite box (e.g. with Dirichlet boundary conditions) - whereas the latter divergence 
5_(r, T) is restricted to the scenario of a black hole, but it is not affected by a finite box. E.g., 
the expectation value of the number "operator" N-p in any KMS state (with a non- vanishing 
temperature) contains the complete divergence (N r ) T ~ 5 + (r,T) + 5_(r,T). One important 
example is the Israel-Hartle-Hawking |9|, state, the KMS state corresponding to the Hawking 
temperature T = k/(2tt). For large radial distances to the black hole the (renormalised) 
expectation value of the energy-momentum tensor evaluated in this state approaches a constant 
value (proportional to T 4 ). In contrast, for the Unruh || state - the state describing the black 
hole evaporation - the (renormalised) energy density decreases with 1/r 2 for large r. As a 
consequence, the expectation value of the number of particles in this state does not display the 
complete divergence <5 + (r,r). 
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It might be interesting to illustrate the point above with the aid of the Bogoliubov coefficients: 
If we consider the spherically symmetric collapse of a star to a black hole the metric outside the 
initial radius of the star does not change (Birkhoff theorem). Ergo the behaviour of the modes 
at very large radial distances r is not affected by the collapse. Accordingly, this region does 
not contribute to the /^^/-coefficients and generates a, S(uj — u/)-term for the a^z-coefficients, 



cf. also P|- Recalling the unitarity relation for the Bogoliubov coefficients in Eq. fl8l|) we 
arrive at the conclusion that exactly this term generates the (5+(T, T)-part of the infinite vol- 
ume divergence. Following Hawking |2] we assume that - for large initial frequencies u - the 
Bogoliubov coefficients are related via Eq. (|76|) in analogy to Section £y| Proceeding in the 
same way as in that Section we observe that the (^-integration of the absolute value squared of 
the /^^/-coefficients is UV-divergent again. But in contrast to the Unruh effect this divergence 
does not contain 5 + (T, T), but only <5_(T, T) owing to the unitarity relation ([H]). Consequently 
the Minkowski vacuum is a KMS state with respect to the Rindler observer, but it does not 
transform into a KMS state during the collapse to a black hole (if we assume the space-time to 
be asymptotically flat and therefore unbounded). Hawking derived the relation (|76f ) only for 
the finally outgoing particles. But even if this relation would hold for both, the (finally) ingoing 
and outgoing particles, the state would still contain less particles than the corresponding KMS 
state. 

As it became evident in the previous considerations, the vicinity of the horizon of a black hole 



displays many similarities to the scenario of the Unruh effect in Section fO|. Indeed, with 
X = K 2 p 2 the black hole metric approaches the Rindler metric in that region 

ds 2 = (n 2 p 2 dt 2 - dp 2 - R 2 dn 2 ) (1 + 0( X )) , (98) 

together with the angular part dVt 2 . This observation motivates an argumentation analogue to 
that at the end of Sec. fO| , cf. also Indeed, it is possible to prove [[UJ that for a state 
fulfilling the Hadamard requirement (among other conditions, see [fT6| ) throughout the complete 
space-time (and in particular at the horizon) the asymptotic expectation values correspond to 
the Hawking temperature. The ground state of the quantum field (the Boulware state) as well 
as every KMS state (with an arbitrary temperature) obey the Hadamard singularity structure 
away from the horizon, see p3| . But only that KMS state that corresponds to the Hawking 
temperature T = k/(2tt) - i.e. the Israel-Hartle-Hawking state - matches the Hadamard 
condition at the horizon. (The same holds true for the Unruh state.) It can be shown that the 
Hadamard condition is conserved during the dynamics of a C°° space-time. Accordingly, if the 
collapse of a star to a black hole can be described by a C°°-metric, the consideration above can 
be used to deduce the Hawking effect. (The Minkowski vacuum is of course also a Hadamard 
state.) However, dropping the assumption of a C°° space-time the situation becomes less clear. 
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If we compare the outcome of this Section with the Minkowski example, we arrive at the 
conclusion that the formation of the horizon causes a bifurcation in a double sense: 
The total Hamiltonian of the field H splits up into two commuting parts H < and H > which 
account for two independent (physical complete) regions r < R and r > R, respectively. 
Before the horizon has been formed there exists only one generalised eigenfunction for every 
given frequency and fixed angular behaviour. This property forbids the definition of ingoing 
and outgoing particles (see Section |4.1|) . After the horizon has been formed the spectrum is 
twice degenerated and the definition of ingoing and outgoing particles becomes possible. 



5.2 Black hole interior 

Our previous investigations focused on the exterior of the black hole. As indicated before we 
shall now take the interior region into account. Inside the (non-extreme) black hole it yields 
h{r) < and therefore g tt < 0, g rr > 0, g&$ < and g w < 0. As a consequence the signature 
of the metric is degenerated and thus the particle and vacuum definition proposed in Section 
|3] does not apply. However, it is still possible to obtain a self-adjoint /C-operator governing the 
dynamics of the system. But for this purpose some modifications are necessary with the result 
that /C is not given by T>'T> and therefore not non-negative. As we shall see later the negative 
parts of K correspond to unstable fields modes. 

At first the scalar product of the interior region {-|-}^ has to be defined with |g 00 | instead of 
g 00 in Eq. (|14D in order to obtain a positive-definite bilinear form. For reasons of simplicity we 
restrict our further considerations to the Schwarzschild metric h(r) = 1 — R/r and start with 
the functions 

C , oo (G<) = lin {C °°(0 <r<R)(g) C°°(§ 2 )} , (99) 

where §2 denotes the 2-sphere of $ and (p. Again the Hilbert space L 2 (G<, |g 00 |) is given by 
the completion of all these functions with respect to the (redefined) scalar product. 
The degenerated signature permits the definition of a scalar product containing g l K Accord- 
ingly, the subsequent steps in Sec. |3| cannot be adopted here. In particular we cannot introduce 
an operator V such that the self-adjoint /C-operator is represented by the absolute value squared 
of T>. Instead we may define an operator /Co via 

K„ : C„°°(G<) c L 2 (G<, | s °°|) - -MG<,|9 M |) 

-f - -£|v£- kv> **- (100) 

The second term at the r.h.s. of the above expression for /Co generates the negative parts of this 
operator. These negative parts originate from the angular derivatives and cannot be obtained 
in a purely radial symmetric consideration. 



30 



Obviously /Co is Hermitian with respect to the scalar product containing the weight \l/h\ (and 
sf^g = r 2 sin#). In addition - since Def(/Co) = Cq°(G < ) is dense in the underlying Hilbert 
space L 2 (G < , \g 00 \) - it is densely defined and therefore symmetric. 

Now we can make use of the following theorem (see e.g. |j3~5|): Every symmetric and C-real 
operator acting on a complex Hilbert space possesses (at least one) self-adjoint extension(s). 
As a result we will always find a self-adjoint operator /C (as an appropriate extension of /Co) 
governing the dynamics of the field. In terms of /C the Lagrange function in Eq. ( |T3"D for the 
interior domain assumes the simple form 

L< = -i i{$|/C|$}< . (101) 

Note that in contrast to Eq. (P0|) the global sign has changed. However, this global sign does 
not affect the equation of motion, but - as it will become evident later - the negative parts of 
the /C-operator do so. 

Since the self-adjoint /C-operator represents an extension of the original operator /Co these two 
operators have to coincide on the subspace Cq°(G < ). Accordingly, it is possible to construct 
test functions w(r, <p) = w(r) D^ m ($, (p) G Cq°(G < ) generating negative expectation values of 
the /C-operator via 

R 

{w\JC\w} = J dr (\h\r 2 \d r w\ 2 -£{£+ 1) |w| 2 ) < 0. (102) 
o 

If we choose the angular quantum number i very large the expectation value {w\ /C \w} 
equals negative numbers of arbitrarily large absolute values, even for normalised test functions 
= 1. Hence the spectrum of /C is unbounded from below. (Of course it is also unbounded 
from above.) Diagonalising the Hamiltonian by means of a quasi-unitary transformation V in 
analogy to Sec. || yields 

H< = ~\^ (P r 2 + Ar^) • (103) 
r 

The interior Hamiltonian is still self-adjoint (by Stone's theorem) - but it is not bounded from 
above and below. Ergo it does not possess a ground state, and a definition of particles as 
excitations over the ground state is impossible. 

As mentioned before, the global sign does not affect the equations of motion, but the occurring 
negative eigenvalues Ar do so: The modes Y corresponding to negative eigenvalues Ar obey the 
following equations of motion 

^Qr = |A r |Qr. (104) 
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Their solutions exp(±y|Ar[i) display a highly (linear) unstable behaviour. 
This instability cannot be avoided by introducing an indefinite metric of the Fock space {13 
if we assume the black hole to be formed by a collapse because in this case the Fock space is 
initially well-defined and obeys a positive definite metric > 0. 

One might suspect that the initial conditions are just in such a way that the exponentially 
increasing solutions do not occur, cf. also [Q. Employing an analogue from classical mechanics 
this situation corresponds to a point-mass moving on the top of a parabolic hill which just comes 
to rest at the zenith of the parabola. However, within quantum theory no (regular) stationary 
state exists in such a scenario (Heisenberg uncertainty relation). Even if the expectation value 
of the amplitude Qr vanishes for all times, its variance increases exponentially (for late times). 
Since the unstable behaviour described above accounts for the time-evolution of the (global) 
modes T it describes a global instability which should not be confused with the concept of local 
stability usually associated with the Hadamard condition, cf. |L3| and |23 



It should be mentioned here that potential terms (which we have omitted in Sec. 0) may also 
give raise to negative parts of the /C-operator. E.g., if the assumptions in Section |T2] fail and 
the scalar curvature 1Z assumes negative values over a large enough volume the operators B 
and /C are not non-negative. However, in this situation the /C-operator is still bounded from 
below (if 1Z does not diverge). Hence only modes up to a certain quantum number are unstable. 
These modes are strongly correlated to the global structure of the space-time. Special care is 
required concerning the interpretation of the instability caused by mass terms. Mass terms that 
are generated by the Higgs mechanism occur in the effective Lagrangian for low excitations and 
cannot be extrapolated to large amplitudes. Restricting ourselves to the massless and minimally 
coupled scalar field (as a model for the photon field) only the instability due to the angular 
derivatives remains where all these objections do not apply. 

In order to interpret the instability it might be interesting to investigate the corresponding 
proper or generalised eigenf unctions. Near the horizon (inwards), the modes behave as 



fr ~ exp (-n VX^j ~ (2k[R - r ]y^/^) . (105) 

Depending on the behaviour at the origin r = one might expect the existence of proper 
eigenfunctions fr at some points of the negative part of the spectrum. 

However, even if no proper and (pointwise/locally defined) generalised eigenfunctions exist, one 



may still construct suitable distributions fr with analogous properties P2| : Considering the 
Schwartz/Sobolev space <5i(er, V) C L 2 (cr,V) of all continuous functions over the spectrum a 
of the /C-operator we may define a Dirac (^-distribution as a linear functional over this space. 
This distribution 5r = S(X, Ar) is then defined within the dual space ©_i(er, V). It represents a 
generalised eigendistribution of the diagonalised /C-operator V/CV^ 6r = M. 5r = \r$r- Hence 
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its spatial representation /p = V^5r exists at least as a distribution over V^&i(a, V) (which 
is dense in L 2 (G < , \g 00 \)) and describes an eigendistribution of /C. The construction described 
above generates non-vanishing eigendistributions fr for all non-singular points (of the spectral 
measure) T of the spectrum a. Since every open interval of a contains non-singular points we 
can always find an appropriate mode T where fr exists. 

Using these eigendistributions fr{f) we can construct solutions of the Klein- Fock-Gordon equa- 
tion of the form 

F T (r,t) = exp (±y/\X^\t) / r (r) , (106) 

if we choose a mode T from the negative part of the spectrum. As a consequence, the equation 
of motion does not only possess unstable solutions - even the degree of the instability v|Ar| 
can be arbitrarily large. In a vivid description one may speak about an explosion interiorly. 
It should be mentioned here that a (partial) negative Hamiltonian, i.e. a (partial) negative 
generator of the time-evolution, is not sufficient for the prediction of an instability. As a 
counter-example we may consider a 1+1 dimensional black hole with ds 2 = halt 2 — dr 2 /h. In 
this situation there are no angular terms and thus the interior as well as the exterior /C-operator 
are both non-negative. Consequently the equation of motion is completely stable. Of course, 
the interior Hamiltonian H < displays a global minus sign, but this does not affect the equation 
of motion 



H = H> + H < = \f J {PI > + ^,>0f>) -\f (Pt,< + ^r,<Ql<) • (107) 
r,> r,< 

Moreover, although the total Hamiltonian is unbounded from above and below, it splits up 
into two independent parts which are bounded. The existence of a horizon is essential for 
this bifurcation. In a flat space-time the Wightman [^l] axioms (spectral condition) demand a 
non-negative generator for stability. 

The Schwarzschild metric ds 2 = hdt 2 — dr 2 /h or ds 2 = hdt 2 — dr 2 /h — r 2 dQ 2 possesses a 
unique analytic continuation to values of r beyond the horizon r < R. In contrast the analytic 
continuation of the Rindler metric to negative values of p does not lead to a degenerated 
signature and complex values of p and/or r do not describe a physical sheet of the space-time. 
As a consequence one observes no instability in the Rindler metric - i.e. the scenario of the 
Unruh effect. 

The notion of the unstable behaviour obtained above refers to the time t measured by an (out- 
side) observer at a fixed spatial distance to the black hole. One might argue that this time 
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coordinate is not capable for describing effects inside the black hole due to the coordinate singu- 
larity at r = R. However, the instability obtained above is not restricted to the Schwarzschild 
coordinates - it occurs in other coordinate systems as well: By virtue of the transformation 

dt -+ dt± dr . (108) 

1-R/r K J 



the metric of the black hole can be cast into the Painleve-Gullstrand-Lemaitre ESL HBI, 071 form 



ds 2 = (l - dt 2 ±2\jjdtdr- dr 2 - r 2 dVl 2 . (109) 

This metric is regular everywhere except at the singularity at r = 0. The transformation of the 
unstable solutions in Eq. ( |106| ) into this coordinate system via Eq. ( |108| ), i.e. t — > t ± S(r), 
merely results in a simple r-dependent factor 

F T (r,t) = exp (v^AH t) /r(r) exp (±^A^ H(r)) , (110) 

while the unstable behaviour persists. The same holds true for the Eddington-Finkelstein 
coordinates with v — t + r* 

ds 2 = (l- j^j dv 2 -2dvdr-r 2 dQ 2 . (Ill) 

Within these coordinates ingoing light rays are simply governed by v = const. Both coordinate 
systems lead to a stationary (but not static) metric, i.e. the evolution parameter still coincides 
with a Killing vector. (This is not the case for the Kruskal metric.) In summary the instability 
of the field equation inside the black hole turns out to be a quite general phenomenon. 

5.3 Back-reaction 



The Eddington-Finkelstein metric in Eq. ( |1 1 1|) allows for a demonstrative visualisation of the 



unstable behaviour: If one emits radially ingoing light pulses in uniform intervals these beams 
are labelled by equidistant values of v. According to the results of the previous Section the 
amplitude of the field $ inside the black hole increases exponentially with rising numbers v of 
the light rays. Hence we may draw the conclusion that the instability is not just an artifact 
caused by an inappropriate description but a physical effect. 

Nevertheless, for an eternal black hole the outside observer is completely causally separated 
from the region of the instability. Hence the interpretation of the unstable behaviour is not 
obvious in that case. But if one considers the possibility of the decay of the black hole (no 
matter whether via evaporation or explosion) and assumes that this decay can be described 
using one of the coordinates above the unstable behaviour should be relevant. (Of course, the 
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assumption of an eternal black hole automatically excludes some of the scenarios where the 
instability may become relevant.) 

In order to investigate the consequences of the instability one has to deal with the back-reaction 
problem. Within all of our previous considerations the quantum field was regarded as a test 
field, i.e. it did not influence the given (externally prescribed) space-time. It is known from 
classical field theory (see e.g. |48| and references therein) that the formation of the horizon 
and the singularity may well be affected by the scalar field $. (For quantum field theory one 
expects that the back-reaction will become important at the Planck scale.) However, Ref. 
[48fl deals with radially symmetric fields only. For that reason the unstable behaviour was not 



obtained there. The correct implementation of the back-reaction of a quantum field has to 
be determined by an underlying theory unifying gravitational and quantum effects. Since we 
have no well-established solution to this problem, we may only speculate about the impact of 
the quantum field on the metric based on physical reasonable arguments. There are several 
possible consequences: 



• The explosion of the complete black hole 

The unstable field modes evolve as exp ( a/ | Ar | [t — r*]). Hence they "reach" after a finite 
period of time the Planck scale vicinity of the horizon, where the classical treatment 
of the gravitation is expected to break down. In that case one might imagine that the 
"wave front" destroys the horizon and thus the complete black hole. (Such an event 
might perhaps be regarded as a toy candidate for the big bang.) In view of arguments 
concerning the time-reversal symmetry there is no obvious reason why the explosion of 
the complete black hole should be impossible. 

As long as there is some matter falling into the black hole its horizon increases. Depend- 
ing on the particular dynamics of the metric this may prevent the "wave front" from 
"reaching" the vicinity of the horizon. But for a static black hole there is no way to avert 
the impact. 

One should be aware that most of the theorems of classical general relativity - e.g. the 
black hole analogues of the laws of thermodynamics - are based on appropriate energy 
conditions, cf. [g3[]. But incorporating the expectation value of the energy- momentum 
tensor of the quantum field these energy conditions do not hold in general. In some cases 
one may employ averaged energy conditions instead, but even the validity of an averaged 
condition is by no means obvious in view of the unstable solutions of the field equation. 



• The prevention of the singularity at r = 
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One might expect that the impact of the instability is at the origin r = much stronger 
than at the horizon r = R. In fact, also those theorems of general relativity that predict 
a space-time singularity after a gravitational collapse are based on energy conditions. 
Accordingly, taking the back-reaction of the quantum field into account, the formation of 
the singularity may perhaps be avoided. Instead one might imagine some kind of quasi- 
oscillations: Impelled by the (exponentially large) amplitudes of the quantum field, the 
matter around the origin blows up, absorbs the excitations of the field, collapses (while 
the field repeatedly evolves exponentially), and eventually blows up again. 

• The field does not affect the metric 

This possibility cannot be excluded within the framework of quantum field theory in 
given (external) space-times. However, the situation of a completely static black hole 
(neglecting the Hawking effect, which is very small for macroscopic black holes) seems to 
be rather strange. In that case the amplitude of the field exceeds the Planck scale after a 
finite period of time (measured by an outside observer). Hence one would expect drastic 
modifications of the space-time. 



5.4 Sonic analogue of black holes 



In 1980 Unruh |L2| discovered a very interesting model for the kinematics of fields in curved 
space-times. He considered the propagation of sound waves in flowing fluids where the effective 
equation of motion assumes the same form as the Klein-Fock-Gordon equation in curved space- 
times. The effective metric depends on the particular flow profile. Many investigations have 
been devoted to this topic during the last years, see e.g. |2(J, the recent work J2"2"| , and references 
therein. 

Before discussing the consequences of the results of the previous section within this scenario 
we shall repeat the basic ideas: The flow of a fluid can be described by its local velocity field 
v, its density g, and the pressure p. The dynamics of the fluid is governed by the non- linear 
Euler equation 

v + (vV)v + -f = f ext , (112) 
g 

if we neglect the viscosity, and the equation of continuity 

g + V(gv) = 0. (113) 

For reasons of simplicity we restrict our further considerations to a constant speed of sound 
c s . This implies the very simple relation between the density and the pressure p = c^g. If we 
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assume an irrotational flow V x v = 0, we may introduce a generating scalar field v = V$. 
Now we linearise the non-linear system of the two equations above around a fixed background 
solution via 

$ = $ + £$i + O(e 2 ), 
v = v + £Vi + 0(e 2 ) , 
p = p + epi + 0{e 2 ) , 

g = Q + s 6l + O{e 2 ). (114) 
This enables us to consider the propagation of small perturbations - i.e. sound waves - within 



a given flow profile. It turns out JT^| that the potential $x of the fluctuations satisfies the 



Klein-Fock-Gordon equation with the effective (acoustic) metric 

^ = ^( C '-"» V °). (115) 



Ergo sound waves in flowing fluids share a lot of interesting features with fields in curved space- 
times. E.g., the surface of transition from subsonic to supersonic flow represents the acoustic 
analogue of a horizon. For a stationary and radially symmetric flow this surface possesses even 
the properties of an event and an apparent horizon. (Unfortunately this scenario exhibits the 
problem of fluid conservation at r = which has to be evaded in some way.) 



Selecting a particular velocity profile v = ±.r^jR/r z it is possible [^(J to simulate a space-time 



which obeys - up to a conformal factor r~ 3//2 - the Painleve-Gullstrand-Lemaitre |45| , [i"6| , f47 
metric in Eq. ( |109| ). According to the results of the previous section the Klein-Fock-Gordon 
equation possesses unstable solutions inside the black hole. Consequently, also the sound waves 
within the supersonic region obey an instability. The conformal factor mentioned above and 
the coordinate transformation in Eq. ( |108| ) do not alter this conclusion - see the remarks in 
the previous Section. 

In contrast to the "real" black hole, where the consequences of the instability are not a priori 
clear (back-reaction problem), there is no possibility to avoid the instability for the acoustic 
black hole models since in that case t denotes the appropriate time also for an inside observer 
and the sound waves affect the fluid directly. 

In the theory of fluid dynamics, such an instability is a well-known indicator for the breakdown 
of the laminar (irrotational) flow, see e.g. ||49|| . I.e., that flow does not represent a stable 
fixed point of the non-linear equation of motion. Accordingly, any small disturbance will grow 
up exponentially until the non-linear regime has been reached. (It should be mentioned here 



that the unstable behaviour obtained above is not a downstream instability, cf. [49 



since 



the perturbation increases exponentially also at a fixed radius r.) In order to investigate 
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the behaviour of the flow after leaving the unstable fixed point - e.g. pattern formation or 
turbulence - one has to consider the non-linear region. For quantum fields in curved space- 
times one expects to reach the non-linear regime at the Planck scale where the back-reaction 
strongly contributes. 

Recalling the outcome of the previous Section the unstable behaviour of the equation of motion 
results from the angular derivatives of the /C-operator. Ergo we may draw the conclusion that 
the quantum field inside the black hole as well as the supersonically flowing fluid favour a 
spontaneous breaking of the radial symmetry, similar to the formation of a vortex in the drain 
of a basin. 

6 Conclusions 
6.1 Summary 

For a minimally coupled, massless and neutral scalar quantum field $ propagating in an 
arbitrary physical complete and causal space-time M that possesses a static metric of non- 
degenerated signature it is possible to perform a particle definition via diagonalisation of the 
Hamiltonian. 

Application of this method to the 3+1 dimensional Minkowski space-time yields the nonex- 
istence of radial ingoing and outgoing particles. For the 1+1 dimensional Rindler metric we 
exactly recover the well-known Unruh effect. 

If we employ the same formalism in order to investigate a black hole the associated space-time 
splits up into two independent domains, inside and outside the horizon, respectively. Within 
the presented approach a particle definition can be accomplished for the exterior region only. 
The quantum field inside the black hole possesses a highly unstable behaviour. The correspond- 
ing Hamiltonian is unbounded from above and below. Accordingly, it is not possible to define 
a vacuum as its ground state and particles as excitations over this state. 

This instability is not a remnant of an inappropriate description but a physical effect. Due to 
our lack of understanding the unification of quantum theory and gravity the consequences of 
this effect are not altogether clear (back-reaction problem). In view of the sonic analogues of 
black holes - where the unstable solutions go along with the breakdown of the laminar flow 
- one might expect that the instability indicates (at least) the breakdown of the treatment of 
quantum fields in given (externally prescribed) space-times. 
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6.2 Discussion 



In order to elucidate the outcome of the formalism presented in this article it might be inter- 
esting to discuss the main statements together with their relations to other approaches: 
As we have observed in Section |4.1| , the particle definition via diagonalisation of the Hamilto- 
nian (equivalent to the energy) does not allow for the introduction of ingoing and/or outgoing 
particles in the Minkowski space-time. The same holds true for more general regular space- 
times. As a consequence, the vacuum coinciding with the ground state cannot be defined as 
that state that is annihilated by the "operators" A? corresponding to purely ingoing (and/or 
outgoing) components exp(— iujv)/r (and/or exp(— iuu)/r), with v = t + r and u = t — r. 
Instead the ground state gets annihilated by "operators" (strictly speaking, operator- valued 
distributions) corresponding to standing waves, i.e. superpositions of ingoing and outgoing 
components with equal weights. Ergo, considering the collapse of a star to a black hole the 
initial ground state cannot be uniquely and consistently defined by the requirement "no ingo- 
ing/incoming particles/radiation". Ref. || states explicitly: Note that we have not defined the 
vacuum by minimizing some positive-definite-operator expectation value (e.g. the Hamiltonian), 
but we have defined the vacuum as the state with no incoming particles. In order to investi- 
gate the relationship of the state defined in this way and the initial ground state additional 
considerations are necessary. 

In contrast to the Minkowski case the ground state of the quantum field in the exterior black 
hole space-time - the Boulware state - has to be defined via demanding that the action of the 
annihilators for both, the ingoing and outgoing modes, yields zero: V^ m : A^ m \ ^b) = 0. 
This fact illustrates the bifurcation caused by the formation of the horizon. 
However, if we assume the black hole to be enclosed by a large sphere with e.g. Dirichlet 
boundary conditions then the definition of ingoing or outgoing particles is impossible again. 
This observation demonstrates manifestly that the particle interpretation is a global concept - 
it may be influenced by objects (e.g. the sphere) at arbitrarily large distances. 
As another difference between the black hole and the Minkowski situation we may recall the 
fact that the /C-operator of the black hole possesses - even in the presence of a finite sphere - a 
continuous spectrum. Due to the additional effective infinity at the horizon the infinite volume 
divergence of the black hole space-time cannot be regularised by enclosing it by a finite box. 
(This regularisation applies only to space-time without any horizon.) 

There are two main interpretations of the Hawking effect: The first view considers the particles 
to be produced by the dynamics of the space-time during the collapse while within the second 
view the radiation is created in a steady rate after the collapse. The observations in Section [5], 
i.e. the splitting of the total Hamiltonian into two independent parts and the diagonalisation of 
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the exterior part by a suitable particle definition (where the number of particles is conserved), 
supports the former interpretation. 

The Hawking effect may be regarded as the verification of the extension of the laws of ther- 
modynamics to objects like black holes. This effect allows us to assign a temperature to the 
black hole via T = 1/(8M) for the Schwarzschild black hole with h — 1 — 2M/r. As a result 
the associated heat capacity of the black hole turns out to be negative: If the mass/energy 
increases the temperature decreases. The classical laws of thermodynamics predict that an 
object obeying a negative heat capacity will be unstable. Accordingly the instability of the 
black hole interior as observed in Section [|]2] might also be regarded as a verification of the 
application of thermodynamics to black holes. 

The consequences of the unstable behaviour of the Klein-Fock-Gordon equation in the interior 
of the black hole cannot be deduced rigorously within the framework of quantum fields in 
(externally prescribed) space-times. The evaluation of the impact of this instability demands the 
knowledge of the back-reaction which has to be determined by a unifying theory. Nevertheless, 
if the underlying theory possesses an evolution parameter corresponding to the Schwarzschild 
time t (or one of the other coordinates discussed in Sec. |5.2|) and contains the treatment of 
quantum fields and external metrics in some limiting case, then one would expect that the 
representation of a black hole also obeys the linearly unstable behaviour. (This would be 
consistent with the frequently adopted interpretation that black holes are highly excited states 
of the unifying theory.) 

For the situation of the acoustic black hole the interpretation of the unstable behaviour is more 
obvious. Without any mechanism preserving (enforcing) the radial symmetry (e.g. effects of 
super-fluids) it is probably impossible to realise the sonic analogue of a black hole experimen- 
tally. 

6.3 Outline 

The particle definition presented in this article is restricted to static space-times. This in- 
cludes the Schwarzschild and the Reissner metric, but not the Kerr space-time describing a 
rotating black hole. Accordingly, further investigations should be devoted to the extension of 
the previous results to stationary metrics. (Without any Killing vector generating the time- 
translation symmetry it is probably impossible to perform a unique and physical reasonable 
particle definition.) 

Another important extension of the provided formalism is given by the incorporation of the 
electromagnetic field 

C = - A F^. (116) 
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The Maxwell theory possesses primary and secondary constraints [32]. These gauge-problems 
have to be solved before the quantisation and the particle definition becomes possible. One 
way to accomplish this - which seems to be suitable to the canonical approach - is the method 



of separation of variables 0]. Nevertheless there is no obvious reason why the main results of 



this article should not persist. The equation of motion of the electromagnetic field is given by 

V, F>» = -L d M <T g va d P A a ) " 77= ^ (V=9 9 W 9 va d a A p ) = . (117) 

For a very rough estimate one may drop the second term, which is related to the longitudinal 
degrees of freedom. The remaining equation possesses unstable interior solutions similar to the 
scalar field scenario. 
The investigation of the Dirac field 

C^^-fDn-mj* (118) 

around charged black holes creates some new kind of problems, see e.g. [^TJ. Similar to the 
Schwinger mechanism in the semi-classical description a tunnelling process is possible. This 
tunnelling probability gives raise to the question of whether a stable vacuum in the quantum 
field theoretical treatment exists. 

Having obtained a linear instability of the linear equations of motion one may ask whether the 
unstable behaviour persists for non-linear equations of motion including interaction terms, for 
instance $ 4 . One might suspect that the non-linear terms generate new stable fixed points 
of the equation of motion - i.e. a non-perturbative stabilisation of the black hole. However, 
in this case the amplitude of the field has to be located at some fixed scale while the (linear) 
instability exists for arbitrary large scales |Ar|. This might be an argument for the dominance 
of the unstable linear contribution in this region. In order to elucidate this point it is necessary 
to consider the scale behaviour of the interacting theory. 

This article considers the propagation of quantised fields in a given (i.e. externally prescribed) 
space-time. To examine how the quantum fields influence the metric one has to deal with the 
back-reaction problem. Within the canonical (operator) quantisation one usually employs the 
renormalised expectation value of the energy-momentum tensor as the source of the Einstein 
equations |T(J , and within the path-integral approach one may integrate out the quantum field 
in order to obtain an effective action (accounting for the degrees of freedom associated with the 
dynamics of the space-time). However, a complete solution to this question probably requires 
the knowledge of the unification of general relativity and quantum field theory. 



41 



Acknowledgement 

The author is indebted to A. Calogeracos, K. Fredenhagen, I. B. Khriplovich, G. Plunien, G. 
Soff, and R. Verch for fruitful conversations and helpful criticism. Discussions with R. Picard 
concerning questions of functional analysis are also gratefully acknowledged. This work was 
partially supported by BMBF, DFG and GSI. 

References 

[1] S. A. Fulling, Phys. Rev. D 7, 2850 (1973). 

[2] S. W. Hawking, Nature 248, 30 (1974); Comm. Math. Phys. 43, 199 (1975). 
[3] P. C. W. Davies, J. Phys. A 8, 609 (1975). 
[4] W. G. Unruh, Phys. Rev. D 14, 870 (1976). 
[5] D. G. Boulware, Phys. Rev. D 11, 1404 (1975). 
[6] R. M. Wald, Comm. Math. Phys. 45, 9 (1975). 

[7] A. Ashtekar and A. Magnon, Proc. R. Soc. Lond. A 346, 375 (1975). 
[8] J. B. Hartle and S. W. Hawking, Phys. Rev. D 13, 2188 (1976). 
[9] W. Israel, Phys. Lett. A 57, 107 (1976). 
[10] R. M. Wald, Comm. Math. Phys. 54, 1 (1977). 

[11] S. A. Fulling, M. Sweeny, and R. M. Wald, Comm. Math. Phys. 63, 257 (1978). 
[12] W. G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

[13] R. Haag, H. Narnhofer, and U. Stein, Comm. Math. Phys. 94, 219 (1984). 

[14] B. S. Kay, Comm. Math. Phys. 62, 55 (1978); Comm. Math. Phys. 100, 57 (1985). 

[15] J. Dimock and B. S. Kay, J. Math. Phys. 27, 2520 (1986); Ann. Phys. 175, 366 (1987). 

[16] K. Fredenhagen and R. Haag, Comm. Math. Phys. 127, 273 (1990). 

[17] B. S. Kay and R. M. Wald, Phys. Rep. 207, 49 (1991). 

[18] R. Brout, S. Massar, R. Parentani, and Ph. Spindel, Phys. Rep. 260, 329 (1995). 

42 



[19] S. Deser and O. Levin, Phys. Rev. D 59, 64004 (1999). 
[20] M. Visser, Class. Quant. Grav. 15, 1767 (1998). 

[21] I. B. Khriplovich, Zh. Eksp. Teor. Fiz. 115, 1539 (1999); Phys. Rep. 320, 37 (1999). 

[22] S. Liberati, S. Sonego, and M. Visser, Class. Quant. Grav. 17, 2903 (2000). 

[23] R. Verch, Comm.. Math. Phys. 160, 507 (1994); H. Sahlmann and R. Verch, e-preprint: 
|math-ph/000202l . 



[24] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space, Cambridge Univ. Press, 
Cambridge, 1982. 

[25] S. A. Fulling, Aspects of Quantum Field theory in Curved Space-Time, Cambridge Univ. 
Press, Cambridge, 1989. 

[26] R. M. Wald, Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics, 
Univ. of Chicago Press, Chicago, 1994. 

[27] Matt Visser, Lorentzian Wormholes: From Einstein to Hawking, AIP Press, Woodbury, 
1995. 

[28] P. K. Townsend, Black Holes, Lecture notes, DAMPT, University of Cambridge, e-preprint: 



gr-qc/9707012 



[29] L. H. Ford, Quantum Field Theory in Curved Spacetime, Lecture notes, to be published in 
the proceedings of the 9th Jorge Andre Swieca Summer School: Particles and Fields, Sao 
Paulo, Brazil, 16-28 Feb 1997; in Campos do Jordao 1997, Particles and fields, 345-388; 
e-preprint: |gr-qc/9707062 . 



[30] M. Bordag, Quantum Field Theory Under the Influence of External Conditions, Teubner, 
Stuttgart, 1996; The Casimir Effect 50 Years Later, World Scientific, Singapore, 1999. 

[31] A. Wipf, Quantum Fields near Black Holes, in Black Holes: Theory and Observation, ets. 
F. Hehl, C. Kiefer, and R. Metzler, Springer, Berlin, 1998; e-preprint: |hep-th/9801025. 



[32] M. Henneaux and C. Teitelboim, Quantization of Gauge Systems, Princeton Univ. Press, 
Princeton, 1992. 

[33] S. W. Hawking and G. F. R. Ellis, The large scale structure of space-time, Cambridge 
Univ. Press, Cambridge, 1973. 



43 



[34] M. Reed and B. Simon, Methods of Modern Mathematical Physics, Academic Press, New 
York, 1972 (Vol. I) - 1978 (Vol. IV). 

[35] J. Weidmann, Linear Operators in Hilbert spaces, Springer, New York, 1980. 

[36] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover publications, 
New York, 1965. 

[37] R. Kubo, Phys. Soc. Japan 12, 570 (1957). 

[38] P. C. Martin and J. Schwinger, Phys. Rev. 115, 1342 (1959). 

[39] Y. Takahashi and H. Umezawa, Collective Phenonema 2, 55 (1975); Int. J. Mod. Phys. B 
10, 1755 (1996); H. Umezawa, Advanced Field Theory, AIP Press, New York, 1993. 

[40] B. S. DcWitt and R. W. Brehme, Ann. Phys. (NY.) 9, 220 (1960); 

J. Hadamard, Lectures on Chauchy's Problem in Linear Partial Differential Equations, 
Yale Univ. Press, New Haven, 1923. 

[41] A. S. Wightman, Phys. Rev. 101, 860 (1956). 

[42] K. Fredenhagen, private communications. 

[43] B. Schroer and J. A. Swieca, Phys. Rev. D 2, 2938 (1970). 

[44] S. A. Fulling, Phys. Rev. D 42, 4248 (1990). 

[45] P. Painleve, C. R. Acad. Sci. (Paris) 173, 677 (1921). 

[46] A. Gullstrand, Arkiv. Mat. Astron. Fys. 16, 1 (1922). 

[47] G. Lemaitre, Ann. Soc. Sci. (Bruxelles) A 53, 51 (1933). 

[48] L. M. Burko, Phys. Rev. D 58, 084013 (1998). 

[49] L. D. Landau and E. M. Lifschitz, Fluid Mechanics, Course of Theoretical Physics Vol. 6, 
Pergamon Press, London, 1959. 

[50] R. Schiitzhold, G. Plunien, and G. Soff, Phys. Rev. A 58, 1783 (1998). 



44 



